Fermi surface reconstruction in strongly correlated Fermi systems as a first order 

phase transition 
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A quantum phase transition in strongly correlated Fermi systems beyond the topological quantum 
critical point is studied within the Fermi liquid approach. The transition occurs between two 
topologically equivalent states, each with three sheets of the Fermi surface. One of these states 
possesses a quasiparticle halo in the quasiparticle momentum distribution n(p), while the other, 
the hole pocket. The transition is found to be of the first order with respect to both the coupling 
constant g and the temperature T. The phase diagram of the system in the vicinity of this transition 
is constructed. 
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Low temperature quantum phase transitions in 
strongly correlated Fermi systems is one of hot topics 
in the condensed matter physics in the last decade. Vari- 
ation of external parameters (pressure, density, magnetic 
field) allows one to shift the transition temperature to 
zero and to obtain the quantum critical point, which is 
associated with divergence of the effective mass M*. In 
the vicinity of this point, low temperature properties of 
the system possess non-Fcrmi-liquid character, i.e. they 
are not described within the conventional Landau theory 
of Fermi liquid. 

At present, experimental information on the quantum 
critical point is available only for three types of strongly 
correlated Fermi systems: i) the inversion layer in MOS- 
FET siHcon transistors in which electrons form a two- 
dimensional (2D) liquidfii^ ii) films of '^He atoms on var- 
ious substrates, iii) metals with heavy fermions.-- 

In nonsuperfiuid homogeneous and isotropic Fermi sys- 
tems, which will be considered in this work, the ratio of 
the effective mass M* to the bare one M reads 



M 



^ gS(p,£) 



(1) 



where e° = p'^/2M — fi, fi is the chemical potential, S 
is the mass operator, and the quasiparticle weight z in a 
single particle state is given by z = [1 — {dT,{p, e) /de)^]^^ 
(index means evaluation of the derivative on the Fermi 
surface). The formula ([Ij allows one to consider two 
scenarios of the quantum critical point. The collective 
scenario is build on a supposition that energy depen- 
dence of the mass operator prevails over its momentum 
dependence due to exchange by critical fluctuations in 
the vicinity of collapse point of the respective collective 
mode, and leads to vanishing of the quasiparticle weight 
z and, hence, to divergence of the effective mass M* just 
at that point f§r^ The topological scenario of the critical 
point assumes z-factor to be finite at that point, however 
dominating momentum dependence of the mass operator 
results in the change of the Fermi surface topologyJ^"— 



The reader can find comparison of these two scenarious 
in Refs. [islfl^ . In this paper, we consider the topological 
scenario of the quantum critical point. 

In this connection, it is worth to note that in accor- 
dance with topological classification^^ of ground states 
of fermionic systems, the basic classes differ by topo- 
logical dimension V of the manifold of nodes of the 
single- particle spectrum e{p) measured from the chem- 
ical potential. Within the same class, we will distin- 
guish states by a number of connected sheets of that 
manifold. All transitions between ground states which 
belong to different topological classes or transitions be- 
tween states with different topology in the same class are 
quantum phase transitions occurring at T = 0. Conven- 
tional nonsuperfiuid homogeneous and isotropic Fermi 
liquid at T = with quasiparticle momentum distribu- 
tion riFL(p) = S{pp — p) belongs to the class for which 
the dimension 2? of the manifold of nodes is less by unity 
than the dimension of the system itself, and this manifold 
is a single connected sheet, i.e. the Fermi surface. 

Violation of the necessary stability condition for Lan- 
dau quasiparticle ground state with the momentum dis- 
tribution riFL(p) serves a signal for its topological recon- 
struction. This stability condition 



5E^2 e{p, [npLb)]) Snpi^ip) dv > 



(2) 



demands positivity of variation of the ground state en- 
ergy E for any admissible variation dnpi^ (p) which satis- 
fies the condition 



2 / SnFi^{p)dv = . 



(3) 



In Eqs. and ([3]), dv denotes an elementary volume of 
the momentum space, and a factor of two means summa- 
tion over two spin projections. The distribution rtFL(p) 
satisfies the necessary condition ([2]) provided the single- 
particle spectrum e{p, [riFhip)]) vanishes only at p = pp- 
In weakly and moderately correlated systems this is true. 
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FIG. 1: Single-particle spectrum and momentum distribu- 
tion of quasiparticles: illustration of topological reconstruc- 
tion scenarios. Panel (a): hole pocket, panel (b): quasiparti- 
cle halo, panel (c): symmetrical three-connected distribution, 
panel (d): fermion condensate. 



However, in process of correlations strengthening with 
the change of external parameters, new nodes of the func- 
tion e{p, [nFL(p)]) can appear and the condition ([2]) is 
then violatedjiSi"— 

Two scenarios of topological reconstruction of the mo- 
mentum distribution npL (p) resulting from this violation 
are known: i) reconstruction within the same topolog- 
ical class, i.e. with no change of the dimension V but 
with change of the number of connected sheets of the 
Fermi surfacej^^"— ii) reconstruction with transition to 
another topological class, i.e. with change of the dimen- 
sion VJ^^— When new nodes of the spectrum appear 
on the same side from the Fermi surface, the distribu- 
tion UflCp) is rearranged to asymmetric three-connected 
momentum distribution which is schematically shown in 
the panels (a) and (b) of Fig. [TJ If the nodes pi, p2, ps 
are arranged in such a way that \pi — P2 1 ^ b2 — | 
the distribution has a form of the hole pocket in the 
filled sphere (we refer to this state as to "H-state), while 
if \pi — P2I ^ \P2 — Pal, one deals with the quasipar- 
ticle halo (P-state). Such reconstruction results in no 
change of the topological dimension V, while the num- 
ber of connected sheets of the Fermi surface appears to 
equal three. If new nodes of the spectrum emerge to both 
sides of the Fermi surface, then together with forma- 
tion of symmetrical three-connected Fermi surface^Sr— 
(see panel (c) of Fig. [IJ, essentially different scenario 
of rearrangement of Landau state is possible, fermion 
condensation^"— shown in the panel (d) of Fig. [TJ In 
this scenario the quasiparticle momentum distribution 
gradually drops within the interval pt < p < pf, and 
the spectrum e{p) identically vanishes within this inter- 
val. Hence the state with fermion condensate turns out 
to belong to the class with the topological dimension of 
manifold of nodes V coinciding with dimension of the 
system. The fermion condensate, revealed and studied 



in details about 20 years ago^^"— acquires a new life 
in these days in a form of topologically protected fiat 
bands, i.e. dispersionless branches of the single-particle 
spectrum with exactly zero energyi ^^i^^ Particularly, pos- 
sibility of existence of surface states with flat band is in- 
tensively discussed)^""— which may be superconducting 
with high transition temperature.^^ 

In this paper, we consider the scenario of the topo- 
logical reconstruction with formation of three-connected 
Fermi surface. We will show that in a topologically rear- 
ranged system, the first order transition between V- and 
■H-states may occur. 

We focus now on the scenario of topological transition 
in which only regions adjacent to the Fermi surface are 
involved. Results of microscopic calculations for 2D liq- 
uid ■^He^^ and for low-density 2D electron gas^^ indicate 
this way of topological reconstruction in these systems. 
For evaluation of the single-particle spectrum e{p) and 
momentum distribution of quasiparticles n{p) we use the 
Fermi-liquid relation^"— 



dejp) 
dp 



in which n{p) — 9{~e{p)) and the quasiparticle interac- 
tion /(p, pi) in the Fermi liquid theory is supposed to be 
known function of momenta. The formula ^ represents 
the nonlinear integro-differential equation for the single- 
particle spectrum e(p). Any numerical algorithm of its 
solution requires use of regularization procedure. This is 
finite temperature that plays a role of a natural physi- 
cal regularizer. Indeed, making use of the Fermi-Dirac 
relation between momentum distribution and spectrum 



nip,T) 



1 



(5) 



allows one to solve Eq. by standard iterative algo- 
rithm. 

We analyze the topological reconstruction in 2D Fermi 
system with a quasiparticle interaction function 



/(P,Pi) 



TT 



1 



M((p-pi)2/9o'-l)'+/3' 



(6) 



with qo ~ 2pp, /3=0.14, which enables one to reproduce 
adequately microscopic calculations^ of single-particle 
spectra of 2D electronic gas at T = on the Fermi-liquid 
side from the quantum critical point. Since the interac- 
tion function depends on the difference p — Pi, Eq. (|4|) 
is integrated to the form 

p2 f 
eb) = ^ " ^ + / -^(P' Pi) "(Pi) '^^i ' C^) 

in which the chemical potential /i is obtained from the 
normalization condition 



(8) 



2 / n{p) dv — p . 
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FIG. 2: Single-particle spectra e(p) and quasiparticle momen- 
tum distributions n(p) evaluated at T = 10~^ in units of 
£% = pp/2M for the model with the interaction ((6]) with 
go = 2pF for different values of the interaction constant g: 
0.176 (panel (a)), 0.180 (panel (b)), 0.190 (panel (c)), 0.218 
(panel (d)) and 0.224 (panel (e)). 
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FIG. 3: Energy per one particle as a function of pi and p2 
evaluated at T = for go = 2pF- 



Single-particle spectrum e(p) and quasiparticle momen- 
tum distribution are evaluated by self-consistent solution 
of Eqs. dS]), ([7]) and ([S|). Rearrangement of the ground 
state of the considered system with increase of the in- 
teraction constant g is shown in Fig. [2j Calculations 
are performed at T = 10~^ep modeling zero temper- 
ature. Irregularity of the spectrum at p > pp distin- 
guished in the panel (a) of this Figure is developed to its 
nonmonotonous behavior which, as g reaches gf, = 0.18, 
results in the bifurcation in the equation e(p) = at 
P = Pb > Pf (see panel (b)) and then, to the topological 
reconstruction with formation of the T'-state (panel (c)). 

The three-connected momentum distribution at zero 
temperature n^lp) — 9{pi — p) — 0{p2 — p) + ^(Ps — p) 
is determined in the functional space by two indepen- 
dent parameters, the third one being obtained from the 



2 2,2 2 

Pi -P2 +P3 =Pf 



(9) 



following from the normalization condition ([5]) . This im- 
plies that the energy functional of the system 

E[n] = 2 J ^'^iP)dv + J fip,pi)n{p)n{pi)dvdvi 

considered within the class of distributions n^ip) is just 
a function of two variables, say, pi and p2- Evaluation 
of the function E{pi,p2) indicates that the momentum 
distribution n{p) obtained by self-consistent solution of 
Eqs. dni), d?]), dll) and shown in the panel (c) of Fig. [5] 
corresponds to the global minimum of this function. 
There is no other local minimum of E{pi,p2) just be- 
yond the topological transition point. However, the situ- 
ation changes with increasing coupling constant, namely, 
a new minimum appears at g ~ 0.21. The relief of the 
function E{pi,p2) at g = 0.218 is shown in the left panel 
of Fig. [3) The deep minimum at p2 > Pf corresponds to 
the ground T'-state, the quasiparticle momentum distri- 
bution and the spectrum of which are shown in the panel 
(d) of Fig. [21 The shallow minimum at p2 < Pf corre- 
sponds to the metastable "H-state which is obtained by 
solving of the set of Eqs. dS]), dT]), dSJ provided the itera- 
tion procedure is started from a state inside the shallow 
well With further increasing of the coupling constant g, 
the "H-state minimum lowers with respect to the 7^-state 
minimum, and both minima equalize at g — 0.22. The 
first-order transition from the 'P-state to the "H-state oc- 
curs at this point, the latter state becomes the ground 
one at g > 0.22. This is demonstrated in the right panel 
of Fig. [3] where the relief E{pi,p2) at g = 0.224 is drawn. 
The deep minimum corresponds to the "H-state shown 
in the panel (c) of Fig. [21 As follows from calculations 
made up to g = 0.26, the P-state keeps on existing as a 
metastable one. 
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FIG. 4: Phase diagram in (go, 3) variables. Capital letters 
denote the state occupying corresponding part of the diagram. 
Two letters mean that the first state is the ground state, while 
the second one is metastable. 

Analysis of metamorphoses of solutions of Eqs. dS]), 
dZj), dSj) with variation of both the coupling constant g 
and the wave vector go allows one to build the phase di- 
agram of the system in these variables which is shown in 
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FIG. 5: Quasiparticle momentum distributions and single- 
particle spectra evaluated in the model . Panel (a) shows 
unstable Landau state, panel (b) demonstrates the "H-state, 
panel (c) — the P-state. Single-particle spectrum for the 
Landau state e(p, [npL]) shifted by the respective difference 
of chemical potentials is shown for comparison in panels (b) 
and (c) by dashed lines. 
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FIG. 6: Energies per one particle in units of e%- evahiated for 
the "H-state (left part of the figure) and the "P-state (right 
part) at three values of the constant /. 



Fig. m At (7o > 1.96pi?, the diagram is arranged sim- 
ilarly to the one considered above for go — 2pi?. Five 
points on it correspond to five solutions shovifn in Fig. (2] 
Arrangement of the phase diagram at go < 1.96^^ is 
different, namely, the three-connected 7^-state emerges 
just beyond the point of the topological transition and 
remains the ground state vsrhile the metastable T'-state 
appears with increasing g. 

It is viforth noting that the first-order transition under 
consideration is not inherent in 2D systems only. Anal- 
ysis for 3D systems shovi^s that an analogous transition 
occurs in 3D as well. 

Why the considered set of equations possesses simulta- 
neously two solutions at fixed parameters, can be under- 
stood with a help of a simplified model with (5-function 
quasiparticle interaction. 3D system is somewhat more 
convenient for this purpose than 2D one since all calcu- 
lations can be done analytically for the 3D case. For the 



(11) 



with go — 2pF, the state with quasiparticle momentum 
distribution nz(j)) at T = has the spectrum 



fe=i 



Single-particle spectra evaluated with use of Eq. 
with go = 2pF, / = 2.0 are displayed in Fig. [5] The 
spectrum e(p; [^ifl]) given by an account of the only term 
in the sum in (jl2p with the boundary momentum pf is 
shown in the panel (a). Due to 0- function on the r.h.s. of 
Eq. p2|) . the spectrum possesses a kink and changes its 
behavior at the point p^^^ — go— Pf = Vf- The necessary 
condition © for stability of the Landau state with the 
quasiparticle distribution nFL(p) is, evidently, violated. 
The spectra e(p; [ns]) shown by solid lines on panels (b) 
and (c) possesses three kinks. If p2 < Pf, the second kink 



is placed at the point p 



(2) 



9o — Vi lying to the right 



of pf- In this case tuning of the chemical potential to 
the condition of conservation of the quasiparticle num- 
ber gives rise self-consistently to the "H-state. Such state 
with the nodes pi — 0.82pF, P2 = 0.96pF, Vz — 1-lpF 
is shown on panel (b) of Fig. [5] together with the spec- 
trum for the Landau state shifted for convenience by the 
difference of the chemical potentials. In case p2 > Pf, 
the point of the second kink is placed to the left of pf, 
and then tuning of the chemical potential gives rise to 
the T'-state. The spectrum of this state with the nodes 
Pi — 0.97pF, P2 = 1.16pF, P3 = l.lSpF is shown on 
panel (c) together with shifted spectrum e{p; [npi]). 

To elucidate which of the two states, H oi V, proves 
to be the ground one, we evaluate the energies of these 
states. Dimensionless energy of the tree-connected state 
per one particle £ = AE/e^N measured from the energy 
of the Landau state is given as follows 



pI 



2/go 

PF 



-p%) 



(5(go;[n3])-5(go;[nFL])). (13) 



Upon not difficult but cumbersome algebra, the structure 
function 



2 f d ] 

S{q; W) = - / n{p + q)n(p) — 
P J (27r 



(14) 



is evaluated analytically. Excluding, say, the variable ps, 
one then arrives at the energy as a function of two vari- 
ables, pi and p2- The condition of its extremum allows 
one to express p2 via pi and reduce the energy to the 



5 



function of a single variable pi. Let qo = 2pF, we intro- 
duce then a new convenient variable x ~ 1 ~ Pi/pf- For 
small values of a;, the energy of the 7^-state equals 



£v{x) 



-a-p (5p 



-S^)0{S-p)x + b'p 



(15) 



where 6-p — f — is an excess of the coupling constant / 
over the critical value /|> = 1 corresponding to the topo- 
logical transition from the Landau state to the P-state, 
a-p = 3/4, b-p — 6. For the H-state, one analogously 
obtains 



(16) 



where S-h ^ f — f^, /-^ = 1 + l/\/2 is a critical value of 
the constant at which the metastable "H-state emerges, 
a-H = 6(^2 - 1), b-H = 9/2 - ^2. Expressions (HI]) and 
pop , formally applicable in the vicinity of the respective 
critical constants, qualitatively describe behavior of the 
energies of V- and H-phases far from and as well. 

As long as S-p < 0, the linear in x term in the en- 
ergy excess of the T'-state over the Landau state equals 
zero and £p{x) = bp x^ > 0. For such values of the 
coupling constant, (5^ is also negative and, hence, both 
terms in £-h{x) are positive. Therefore, the Landau state 
with the quasiparticle distribution n-pi^ip) is the ground 
state. The functions £p{x) and S-h{x) at / < are 
shown by dotted lines in Fig. [HI At (Jp > 0, linear in 
X term with minus sign emerges in £p{x). As a result, 
the P-state wins the contest against the Landau state. 
Thus, the second order topological transitions occurs at 
f = fp, namely, the three-connected P-state with the 
quasiparticle halo appears to be the ground state. The 
energy S-uix) remains a monotonically growing function 
of X as long as 5-h < 0. Both energy curves at the con- 
stant / = 1.6 corresponding to the case fp, < f < 
are shown by dashed lines in Fig. [SI As soon as / ex- 
ceeds /.^, the coefficient near the quadratic term in the 
function £-h{x) changes the sign and the function acquires 
the minimum, i.e. the metastable "H-state appears. When 
with incresing /, this minimum shown by a solid curve 
in the left part of Fig. [6] becomes deeper than the mini- 
mum of the right solid curve £p [x) , the three-connected 
"H-state with the hole pocket, topologically equivalent to 
the T'-state, becomes the ground state. The transition 
between the V- and H-states is first order. 

We return now to the model of 2D electron gas with 
the quasiparticle interaction The height of the en- 
ergy barrier which, as seen in Fig.[3l is of order 10~^e^, 
determines the scale of temperature at which on can ex- 
pect a transition between %- and T'-states with increas- 
ing temperature. Such transition, indeed, occurs. This 
is caused by the fact that the quasiparticle halo is more 
narrow than the hole pocket, and hence, "melts" faster as 
temperature increases, what is demonstrated in panel (a) 
of Fig. [T] As a result, the entropy Sp{T) of the 7^-state 
increases faster with heating than the entropy S-h (T) (see 
panel (b)), while the free energy F-p{T) decreases more 
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FIG. 7: Panel (a): momentum distributions n(p) for the V- 
and the "H-states at T = 2 ■ 10~'^e'^. Entropy per one particle 
S/N (panel (b)), free energy F per one particle F/N in units 
of e%- (panel (c)) and density of single-particle states 11 in 
units of M/n (panel (d)) as functions of temperature T in 
units of e'^. Parameters g — 0.224, qo — 2pF are used in the 
calculations. 



rapidly than F-u{T) (panel (c)). Both minima equalize 
at Ti ~ 1.2 • lO^'^e^^-, and the first order transition occurs 
from the "H-state to the 7^-state. The entropy and the 
density of states undergo a jump at T = Ti (see panels 
(b) and (d)). The considered transition may have rela- 
tion to observed low-temperature anomalies in specific 
heat and magnetic susceptibility of metals with heavy 
fermionsj^i^i 

In conclusion, we analyzed the reconstruction of the 
Fermi surface of the uniform Fermi system with increas- 
ing coupling constant of the quasiparticle interaction and 
found that the topological transition, in which two new 
connected sheets of the Fermi surface appear, is followed 
by the transition between two topologically equivalent 
states. The Fermi surface of both these states consists 
of three connected sheets, but one of these states, the 
T'-state, possesses a structure of the quasiparticle halo, 
while the second one, the "H-state, that of the hole pocket. 
The transition from the T'-state to the 'H-state is of the 
first order with respect to the coupling constant g. As 
the temperature T increases, the inverse first order tran- 
sition from the H-state to the 7^-state occurs due to 
more rapid "melting" of the narrow quasiparticle halo 
and more rapid increase of its entropy with heating than 
increase of the entropy of the hole pocket. 
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